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Abstract. I study the behaviour of the maximum rms fractional amplitude, rmax and the maximum 
coherence, Qmax, of the kilohertz quasi-periodic oscillations (kHz QPOs) in a dozen low-mass X- 
ray binaries. I find that: (i) The maximum rms amplitudes of the lower and the upper kHz QPO, 
r^ax and rjjjaxj respectively, decrease more or less exponentially with increasing luminosity of the 
source; (ii) the maximum coherence of the lower kHz QPO, 2max' first increases and then decreases 
exponentially with luminosity; (iii) the maximum coherence of the upper kHz QPO, 2max' is niore or 
less independent of luminosity; and (iv) rmax and Qmax show the opposite behaviour with hardness 
of the source, consistent with the fact that there is a general anticorrelation between luminosity 
and spectral hardness in these sources. Both r^ax and 2max in the sample of sources, and the 
rms amplitude and coherence of the kHz QPOs in individual sources show a similar behaviour 
with hardness. This similarity argues against the interpretation that the drop of coherence and rms 
amplitude of the lower kHz QPO at high QPO frequencies in individual sources is a signature of the 
innermost stable circular orbit around a neutron star 
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INTRODUCTION 

Kilohertz quasi-periodic oscillations (kHz QPOs) in the X-ray flux of low-mass X-ray 
binaries have drawn much attention since their discovery, about ten years ago. The 
reason for this continued interest is that since they most likely reflect the motion of 
matter very close to the neutron star (or black hole), these QPOs may provide one of 
the few direct ways of measuring effects that are unique to the strong gravitational- 
field regime in these systems. Often two simultaneous kHz QPOs are seen in the power 
density spectra of low-mass X-ray binaries, the lower and the upper kHz QPO according 
to how they appear sorted in frequency. 

Most of the work on QPOs in these years [see 1, for a review] has focused on the 
frequencies of these QPOs, because those frequencies provide insights into the dynamics 
of the system. For instance, different from the Newtonian theory of gravitation, in 
general relativity the effective potential as a function of radial distance to the central 
source has a maximum. A particle in a circular orbit at that radius would be in unstable 
equilibrium; if perturbed, the particle would fall onto the central object. This radius 
defines the innermost stable circular orbit, or ISCO. No stable orbit around the central 
object is possible inside the radius of the ISCO, which in the Schwarzschild case 
(non-rotating central object) is nsco = 6GM/c^ [2]. If (at least) one of the QPOs 



reflects Keplerian motion around the neutron star [3, 4], the ISCO would set a limit 
to the maximum frequency that the QPOs can reach [3]. By measuring that maximum 
frequency it would be possible to infer the mass and constrain the radius of the neutron 
star. 

From the beginning there has been interest on the other two properties of the kHz 
QPOs, their amplitude and coherence [e.g., 5], but systematic studies of those other 
QPO properties only started to take off slightly later [e.g. 6-9]. 

Recently, Barret et al. [10-12] carried out a systematic study of the kHz QPO co- 
herence and rms amplitude in three X-ray binaries, 4U 1636-53, 4U 1608-52, and 4U 
1735^4. They find that in all three sources the coherence and rms amplitude of the 
lower kHz QPO increase slowly with frequency, and after the coherence and rms ampli- 
tude reach their maximum values, they decrease abruptly as the QPO frequency keeps 
on increasing. (The sudden drop is most noticeable in the coherence of the lower kHz 
QPO; in the case of the rms amplitude the drop is less abrupt.) Barret et al. propose 
that this behaviour is due to effects related to the ISCO around the neutron star in these 
systems. 

Triggered by these results, I investigated the dependence of the maximum coherence 
and rms amplitude of both kHz QPOs in a large sample of low-mass X-ray binaries. 

DATA SELECTION AND ANALYSIS 

All the data that I use here were obtained over the last 10 years with the Proportional 
Counting Array (PCA) on board the Rossi X-ray Timing Explorer (RXTE). I collected 
most of these data from the literature. For this I searched, for as many sources as 
possible, all published values of the rms fractional amplitude and the coherence of the 
kHz QPOs. [see 13, for a detailed description of the selection and characteristics of the 
data.] The coherence 2 of a QPO is defined as Q = Vgpo/A, where Vqpo and A are 
the frequency and the full- width at half-maximum of the QPO. The rms amplitude, r, is 
calculated from P, the integral from to oo of the Fourier power under the QPO, and the 
source intensity, 5, as r = 100 x \/P/S; from this definition, r is expressed as a percent 
of the total intensity of the source. 

Here I use the naming convention of [14], in which the lower kHz QPO is called 
Li, and the upper kHz QPO is called L^. The frequency, coherence and amplitude of 
each QPO carry a subscript or superscript £ or u, respectively. I use the expression 
"kilohertz QPO" to refer to features in the power spectra of neutron star systems that 
have frequencies > 150 Hz, and that have not been identified as hectohertz QPO [e.g., 
8]. 

Once I have collected all Qmax and rmax values for each QPO for each source, and the 
frequencies Vi and Vu at which those maximum values occur, I use Figure 1 in [15] to get 
the corresponding source luminosity: Using the QPO frequency as input, I read off the 
luminosity of the source at that frequency from that Figure. For three sources, 4U 1608- 
52, 4U 1636-53, and 4U 1820-30, the maximum rms amplitudes of La occur in states of 
the source in which the QPO frequencies are outside the range of frequencies in Figure 
1 of [15]. In those three cases I either search the literature for flux measurements of the 
source in that state, or I extract spectra from the corresponding RXTE observations, and 




FIGURE 1. Left panel: The maximum coherence (upper panel) and maximum rms amplitude (lower 
panel) of the lower kHz QPO as a function of the source luminosity at the time at which those maximum 
values were reached. The luminosity is in units of the Eddington luminosity for a 1 .QM© neutron star. 
Filled symbols indicate measurements over the full energy band covered by the PCA on board RXTE. 
Open circles indicate measurements over a limited energy band; in these cases (except for Sco X-1), the 
rms amplitudes have been divided by 1.25 [see 13, for details]. The rms amplitude in the case of Sco 
X-1 is not corrected for dead-time, and hence is only a lower limit (indicated with an arrow pointing 
upwards in the lower panel). Circles are Atoll sources; squares are Z sources. Right panel: The maximum 
coherence (upper panel) and maximum rms amplitude (lower panel) of the upper kHz QPO as a function 
of the source luminosity at the time at which those maximum values were reached. The luminosity is in 
units of the Eddington luminosity for a 1 .9M© neutron star. Symbols are the same as in the left panel. 



calculate the luminosities myself, or both. 

The luminosities used here are uncertain for the following reasons: (i) Statistical error 
in the fluxes derived from model fitting. These errors are usually less than 5 — 10%. (ii) 
Accuracy in the determination of the luminosity for a given QPO frequency using the 
results of [15]. The reason for this is the so-called parallel-track phenomenon [e.g., 16], 
which introduces an uncertainty of the order of 20 — 50% in the luminosity [9]. (iii) 
The use of the 2 — 50 keV flux as a measure of the bolometric flux of the source. This 
effect contributes uncertainties of the order of 20 — 25% [see 15]. (iv) Uncertainty in the 
distance to these sources [see 15, for a discussion in the context of these sources]; this 
can yield uncertainties of up to 60% in the luminosity [see, e.g., 17]. In this paper I use 
a fixed error of 25% in the values of L/L^^ as indicative of the error in the luminosity. 
It is clear from the previous discussion that this is a lower limit to the real error in this 
quantity. 

RESULTS 

Figure 1 shows the dependence of the maximum coherence and maximum rms ampli- 
tude of both kHz QPOs as a function of source luminosity (in units of the Eddington 
luminosity for a 1.9Mq neutron star). From this Figure it is apparent that r^^^ and r^^^ 
both decrease more or less exponentially with L/L^dd- Using roughly the same sam- 
ple of sources, [18] had already noticed that the rms amplitude of the upper kHz QPO 
decreases as the luminosity of the source increases. 
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FIGURE 2. The maximum coherence (upper panels) and maximum rms amphtude (lower panels) of 
the lower kHz QPO (left) and the upper kHz QPO (right) as a function of the average source hardness, HI, 
defined as the ratio of the 40 - 80 keV to the 13 - 25 keV count rate, measured with HEAO-1 [25, 26]. 
For an explanation of the symbols see the caption of Figure 1. The rms amplitudes in the case of Sco 
X-1 are not corrected for dead-time, and hence are lower Umits (indicated with a vertical arrow in the 
lower panels). Upper limits to the hardness are indicated with horizontal arrows pointing to the left. I do 
not include GX 5-1 in this Figure because the HEAO-1 measurements of this source suffer contamination 
from a previously unknown hard X-ray source in the field [see explanation in 26]. 



From tlie same Figure, it is also apparent that at low luminosity, Qmax A^'st increases 
with L up to L/ LEdd ~ 0.04, and then it decreases exponentially. On the other hand, Q^ax 
does not show any significant trend with luminosity. 

The gap in Figure 1 at L/Lgdd ~ 0.25 — 0.7 separates the Atoll sources at low L, and 
the Z sources at high L [see 19, for a definition of Atoll and Z sources]. That gap would 
be occupied by the four intermediate-type sources, GX 9+1, GX 9+9, GX 3+1, and GX 
13+1, which so far have not shown any kHz QPOs [20-23]. The upper limit to the rms 
amplitude of the QPO in these sources ranges from 1.6% to 2.6%. Since the range of 
luminosities spanned by these sources [17] is L/L^^ ~ 0.12 — 0.44, these upper limits 
are much lower than would be expected from the interpolation of the trends of r^^^ 
Cax with L/Lfidd in Figure 1. 

From Figure 1 it is apparent that the maximum rms amplitude of both kHz QPOs 
and the maximum coherence of the lower kHz QPO are consistently lower in the Z 
sources than in the Atoll sources (see also below). Since the detection of QPOs in Z 
sources generally require calculating average power spectra over longer time intervals 
than for Atoll sources, this difference could in principle be due to the frequency drift 
of the QPOs, which would artificially reduce their Q values. However, measurements of 
the coherence of the lower kHz QPO in the Atoll source 4U 1608-52 and the Z source 
Sco X-1 on short time intervals (less than a few hundred seconds), when Vf ~ 600 Hz, 
yields Qe = 74.0 ±4.6 for 4U 1608-52, and Qi = 4.2 ±0.4 for Sco X-1. This shows 
that in the Z source Sco X-1, even over very short time intervals, the lower kHz QPO is 
intrinsically significantly broader than in the Atoll source 4U 1608-52. 

From Figure 1 [see also 13] it is also apparent that, except for the case of 4U 0614+09, 
2max ^nd r^^^ are positively correlated with each other. Concerning 4U 0614+09, the 
hardest source in the sample and the one at the lowest luminosity, it is as if in this case 
Gmax were too low for r^^^' or conversely, as if were too high for Q^^. Recent 



work [24] confimis that the decrease of Q^^^ at low luminosity is real. 

Van Paradijs and van der Klis [25] have shown that there is a general correlation 
between the average source luminosity and the average source hardness, HI, defined 
as the ratio of the count rate in the 40 — 80 keV band to the count rate in the 13 — 25 
keV band, measured with HEAO-1 [see 26]. In that respect, as expected, the plots of 
Gmax and Tmax VS. hardncss in Figure 2 show that and both increase with the 
hardness ratio HI, Q^^x increases with HI and then it decreases for 4U 0614+09, the 
hardest source in this sample, and as in the plots as a function of luminosity, Qmax is 
consistent with being constant with HI. I caution the reader that contrary to the L/L^dd 
values that I plot in Figure 1, for this Figure I use average values of the spectral hardness 
measured several years before the kHz QPOs were discovered. Notice that GX 340+0, 
for which the distance, and hence the luminosity, is uncertain, appears in Figure 2 close 
to the other Z sources (open symbols; GX 340+0 is the point at HI = 0.036, just to 
the left of GX 17+2 at HI = 0.037); since HI is a distance-independent parameter, this 
suggests that the distance to GX 340+0 is not too much in error. 

It is interesting to notice in this Figure that there is no gap in the distribution of 
sources as a function of hardness between the Z and Atoll sources. This is opposite 
to what is apparent in the plot of QPO parameters vs. luminosity, where there is a gap 
corresponding to the intermediate-type, the GX-sources (see above). The HI values of 
the GX-sources range from 0.03 to 0.14. 

Both in Figure 1 and 2 there appears to be a gap in the values of Q^^x' ^nd perhaps 
also r^^x and rj^^^ between the Z and Atoll sources. I have just shown that this difference 
is real, and is not an artifact due to the way the QPOs are measured. While this difference 
may indicate a dependence on luminosity (see Fig. 1) or on spectral hardness (see Fig. 2), 
this could also point to a difference between Z and atoll sources. It is worth noting, 
however, that there is still a trend of Q^^x both rms amplitudes within the atoll 
sources in Figure 1. Furthermore, there is a significant trend in the relations between 
Qmax and Tmax within the atoll sources; e.g., the relation between Q^ax ^nd r^jjx is 
8(7 different from a constant [see 13, for details]. All this suggests that the distinction 
between Z and Atoll sources cannot be the (only) explanation for this difference. 

DISCUSSION 

I study the maximum amplitude, rmax, and maximum coherence, 2max, of the kHz QPOs 
as a function of luminosity and hardness for a large sample of low-mass X-ray binaries. 
I show that the maximum coherence of the lower kHz QPO, Q^ax' ^^^^ increases up to 
L ~ 0.04LEdd and then decreases with luminosity, whereas the maximum coherence of 
the upper kHz QPO, 2max' is independent of luminosity. I also find that the maximum 
rms amplitudes of both the lower and the upper kHz QPOs, and r}Jj^, respectively, 
decrease monotonically with luminosity and increase monotonically with the hardness 
of the source. [The dependence of on luminosity and hardness was first reported by 
18]. 

From the above results it follows that for all sources, and rjjj^ are positively 
correlated with each other. Also, for all sources with L>0.04LEdd5 that is all sources in 
this paper except 4U 0614+09, the hardest source in the sample and the one at the lowest 



luminosity, Qj^ax positively correlated both with r^^^ and rjjj^. Qmax independent 
of 2^ax or the maximum rms amplitude of the kHz QPOs. 

In individual sources, both and increase with and then drop rather abruptly 
at the high end of the Vi range; r-a also increases and then drops at high Vu values, and 
2u is more or less constant or increases slightly with Vu [e.g., 7, 9, 27]. In the case of 
4U 1636-53, [11] interpret the sudden drop of the coherence and rms amplitude of L^, 
together with the existence of a frequency above which is not detected, as evidence 
of the innermost stable circular orbit around the neutron star in this system. 

From the results of individual sources and those of the sample of sources that I present 
in this paper, it is apparent that the behaviour of the coherence and rms amplitude of the 
kHz QPOs as a function of the QPO frequency in individual sources is similar to the 
behaviour of the maximum coherence and maximum rms amplitude of the kHz QPOs 
as a function of luminosity in the sample of sources. Since in individual sources there 
is a general relation between QPO frequencies and source intensity, in the sense that 
at higher intensity the QPOs generally appear at higher frequencies (but remember the 
parallel-track effect), this raises the question of whether the same mechanism may be 
behind both behaviours. 

The link between these two behaviours need not be luminosity, but could be the high- 
energy emission (or hardness) in these systems. On one hand, in the sample of sources 
the maximum coherence and maximum rms amplitude of the kHz QPO, except the 
maximum rms amplitude of La, appear to correlate fairly well with spectral hardness 
(see Figure 2), while on the other hand in individual sources the QPO frequencies are 
well correlated with the spectral hardness of the source [16], the index of the power law 
that fits the high-energy part of the X-ray spectrum [28], or S, a parameter that measures 
the position along the track traced out by the source in a colour-colour or colour-intensity 
diagram [called and Sa in the Z and Atoll sources, respectively; see e.g., 6, 29]. From 
this, it follows that in individual sources there should be a relation between QPO rms 
amplitude and coherence on one hand and spectral hardness on the other. 

The comparison between individual sources and the sample of sources suggests that 
the same mechanism is responsible for the drop of coherence and rms amplitude of the 
lower kHz QPO with QPO frequency in individual sources as well as for the drop of 
maximum QPO coherence and maximum QPO rms amplitude with luminosity in the 
sample of sources. Most likely the mechanism is related to the high-energy emission 
in these systems. This does not necessarily mean that the fractional emission at high 
energies (represented by the hardness or X-ray colours) is the root mechanism that drives 
all QPO parameters (QPO frequency, coherence, and rms amplitude). For instance, one 
possibility (there could be many others) is that the (instantaneous) mass accretion rate 
sets the size of the inner radius of the disc [30], which in turn determines the QPO 
frequency as well as the relative contribution of the high-energy part of the spectrum to 
the total luminosity. If the efficiency of the modulation mechanism (related to the rms 
amplitude r) and the lifetime of the oscillations (related to the coherence Q) that produce 
the QPO depended upon the emission from the high-energy part of the source spectrum 
[see 13, for a discussion of possible ways in which this could happen], observationally 
it would appear as if the coherence and rms amplitude of the QPO depended upon the 
QPO frequency, and hence upon the radius in the disc at which the QPO is produced. The 



sudden drop of the coherence and rms amplitude of the QPO at some QPO frequency 
would then appear to be associated to a dynamical peculiarity in the accretion disc, for 
instance the ISCO. Observing the same source repeatedly would not allow to distinguish 
the above scenario from one in which QPO coherence and rms amplitude were actually 
set by QPO frequency or the dynamics in the accretion disc. 

To distinguish between these two options, one would need to observe a sample 
of sources of kHz QPOs for which the mass-accretion rate, and hence the relative 
contribution of the high-energy part of the spectrum to the total emission, was different 
from source to source. In that case, QPO coherence and rms amplitude would drop for 
sources accreting mass at higher rates, even if the frequency of the QPO was more or 
less the same from source to source. Since, as I have shown in this paper, this is the 
general behaviour observed in sources of kHz QPOs, it is reasonable to infer that a 
mechanism similar to the one I described in the previous paragraph is effective in setting 
the coherence and rms amplitude of the kHz QPOs. If this is correct, this also implies 
that the drop of QPO coherence and rms amplitude as a function of QPO frequency in 
individual sources cannot be due to effects of the ISCO. 

Note also that in individual sources not just the rms amplitude of the kHz QPOs, but 
also the rms amplitude of other lower-frequency QPOs decrease with increasing QPO 
frequency. For instance, in four Atoll sources, 4U 1728-34, 4U 1608-52, 4U 0614+09, 
and 4U 1636-53, the rms amplitudes of the "bump", a QPO at ~ 0.1 - 30 Hz, the 
"hump", a QPO at ~ 1 - 40, and the hectohertz QPO at - 100 - 300 Hz, all drop as 
the frequencies of the kHz QPOs increase, in a similar fashion as the amplitude of the 
upper and lower kHz QPOs [e.g., 31, and references therein; see also there a description 
of these other QPOs]. This also argues against the interpretation of the ISCO as the 
cause of the drop of the rms of the kHz QPOs, and indicates that the amplitudes of all 
variability components are set by the same mechanism which, as I suggested, could be 
the same one that governs the high-energy spectral component. 



CONCLUSION 

I study the maximum amplitude, rmax, and maximum coherence, 2max, of the kHz 
QPOs as a function of luminosity and hardness for a dozen low-mass X-ray binaries. 
I find that: (i) The maximum coherence of the lower kHz QPO, Q^^^, first increases 
up to L ~ 0.04LEdd and then decreases with luminosity, (ii) The maximum coherence 
of the upper kHz QPO, QJJ^ax' independent of luminosity, (iii) The maximum rms 
amplitudes of both the lower and the upper kHz QPOs, r^^^ and rj^^x' respectively, 
decrease monotonically with luminosity, (iv) Both and r^ax increase with the source 
hardness, Q^ax ^^^^ increases with hardness and then drops for the hardest source in the 
sample, and Qmax is independent of hardness, (v) The relation between Qmax and rmax 
with luminosity in the sample of sources is similar to the relation between Q and r with 
QPO frequency in individual sources, (vi) The above argues against the interpretation 
that the drop of QPO coherence and QPO rms amplitude at high QPO frequency in 
individual sources is due to effects related to the innermost stable orbit around the 
neutron star in these systems. 
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